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The resistivity around the ferromagnetic transition tem- 
perature in the double exchange model is studied by the 
Schwinger boson approach. The spatial spin correlation re- 
sponsible for scattering of conduction electrons are taken into 
account by adopting the memory function formalism. Al- 
though the correlation shows a peak lower than the transition 
temperature, the resistivity in the ferromagnetic state mono- 
tonically increases with increasing temperature due to a vari- 
ation of the electronic state of the conduction electron. In the 
paramagnetic state, the resistivity is dominated by the short 
range correlation of scattering and is almost independent of 
the temperature. It is attributed to a cancellation between 
the nearest-neighbor spin correlation, the fermion bandwidth, 
and the fermion kinetic energy. This result implies the impor- 
tance of the temperature dependence of the electronic states 
of the conduction electron as well as the localized spin states 
in both ferromagnetic and paramagnetic phases. 



The j-ejueat discovery of colossal magnetoresistance 
(CMR)El'cruu has revived interest in perovskite mangan- 
ites such as Lai-zSr^MnOs. It is widely accepted that 
significant changes in the transport properties, as well 
as CMR, are observed around the transition between the 
ferromagnetic phase and .-the paramagnetic one. More 
than 40 years ago, ZenerH proposed a double-exchange 
(DE) interaction to explain the correlation between elec- 
trical conduction and the ferromagnetism, in which the 
spin of a conduction electron and a localized core spin 
(S) on the same site are strongly coupled by Hund's 
rule. Since the hopping amplitude of the electron to the 
neighboring sites is maximum when the two neighboring 
core spins are parallel, the ferromagnetic metallic state is 
achievecLby-gaining the kinetic energy of the conduction 
electronfrElQ These concepts were settled as the so-called 
DE model and the magnetic and transport properties in 
this model have been investigated intensively and exten- 
sively. 

One of the main interest in this research field is the 
temperature dependence of the electrical resistivity. The 
resistivity in the DE model was studied in a mean-field 
theory by Kubo and Ohatajj and similar results have 
been reproduced by a dynamical mean-field theory by 



Furukawa.EI In these calculations, however, the spatial 
correlation of the core spins is not included properly, al- 
though it is pointed out that it plays a crucial role ii)-4h*| 
electric transport near the Curie temperature (T c )113'til 
The short range spin correlation was only considered in 
Ref. ||, and the spatial correlation was neglected in Ref. ||, 
where the dynamical fluctuation was taken into account. 
Millis et al. discussed a possibility that the behavior of 
the resistivity is greatly modified when the spatial corijfe. 
lation of the core spins is properly taken into account S3 
They showed that the resistivity is still increased below 
T c with decreasing the temperature. Being based on the 
calculated results which disagree with the experimental 
one, they concluded that the additional gradients, such 
as, the Jahn- Teller effect, are necessary to reproduce the 
observed behaviors. 

In this paper, we calculate the temperature depen- 
dence of the resistivity in the DE model by the Schwinger 
boson approach. In order to include the effects of the spa- 
tial correlation of the core spins properly, we adopt the 
memory function formalism which was also used in Ref. 
[l^ . In addition, the temperature dependence of the elec- 
tronic structure is determined self-consistently together 
with that of the core spins. The calculated resistivity 
monotonically decreases with decreasing temperature in 
the ferromagnetic states and does not show a peak be- 
low T c , although the spin correlation has its maximum 
at T < T c . 

The Hamiltonian of DE model in the limit of strong 
Hund's coupling is given by the Schwinger boson repre- 
sentation as follows, 



h.c. 



(1) 



with the local constraint b ia bi a — fjfa — 2S at every 
lattice site i. Here, b ia (cr =|, |) is a boson jand.il is a 
spinlcss fermion operator, Sr =S + (1 — a;)/2,E3o'Ej and 
x is the doping concentration of holes ((///«) = l—x). We 
shall exclusively consider the case of S = | . A transition 
at T c from a ferromagnetic state to a paramagnetic state 
(described as Bose condensation of Schwinger bosons) 
was investigated by using a mean field Hamiltoniantj 
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with a global constraint 

(J2blb ia ) = 2S R , 



(2) 



(3) 



where B and D are given by \Y,A b \<7 b 3°) and (/i/j)> 
respectively, and both are determined self-consistently. 
This mean field treatment, however, leads to an addi- 
tional transition at slightly higher temperature than TL 
(about 1.4T C ) into an artifact state in which B — D — 0X3 
As a result, above T c the fermion bandwidth (W/ = 
l2Bt/Sit) rapidly decreases. This decrease obviously 
causes a misleading diverging increase of the resistivity. 

In the present study, we assume that Eq. (||) itself has 
a suitable form as a mean field Hamiltonian, but in or- 
der to avoid the difficulty mentioned above, the fermion 
bandwidth B is determined as 
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together with 
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in a self-consistent manner. Here, we ignore the Berry's 
phase in the electron hopping and use the fact that 

Y^aa' b io b i°' b \<x' b io ~> 2 ^fl in tlie high-temperatures 
limit. The approximation in Eq. (^) corresponds to the 
expansion with respect to Si ■ Sj. It should be noted 
that the fermion bandwidth obtained in Eq. (4) remains 

finite for T-xx) as expected, while | ^2 a ( b ia b j<r) char- 
acterizing a nearest-neighbor magnetic correlation (de- 
noted by C) vanishes for T— >oo in this model. Therfore, 
the behavior obtained in the above formulas is physically 
reasonable. 

The Curie temperature (T c ) in this model as a function 
of the doping concentration is shown in the inset of Fig. 
[I]. Since T c is determined by Eq. (||), the results are 
essentially independent on the formula of B and are the 
same as those in Ref. |l4]. The transition temperature is 
scaled independent on x by D as T c /t^2.5D. 

The temperature dependencies of the fermion band- 
width normalized by a bare band-width (W = 12t) are 
plotted in Fig. |l|. At the transition temperature, the 
bandwidth does not directly depend on the doping con- 
centration (it depends on x only through Sr). This is 
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FIG. 1. The fermion bandwidth normalized by bare band- 
width (Wf/W) and the fermion kinetic energy normalized by 
the zero-temperature value (Kf /Kf(T = 0)) as a function of 
T/T c for several doping concentrations. The inset shows the 
transition temperature from a ferromagnetic state to a para- 
magnetic state as a function of the doping concentration. 



because both the chemical potential and the condensate 
density of bosons are zero, and Dt/T c , which is inde- 
pendent of x, is only a parameter in Eq. (||).t-3 Further, 
Wf/W approaches 1/V2 in an infinite temperature limit. 
As a result, the behavior of the fermion bandwidth as a 
function of T/T c is almost universally independent of x. 
On the other hand, the behavior of the kinetic energy of 
the fermion (K f = Bt(f}fj)/S R ) as a function of T/T c 
(Fig. |l|) changes depending on x. It is important that 
the bandwidth varies even in the disordered-spin regime 
of T > T c . This behavior agrees well with the result in 



Ref. 16. In fact, the bandwidth at T c is 1.16 times bigger 
than that in the T— >oo limit. 

The resistivity as a function of the temperature is cal- 
culated by the memory function method,E2lL3 where the 
lowest order fluctuation from the mean field can be in- 
cluded automatically. In this lowest-order perturbational 
treatment, a static approximation for Schwinger bosons is 
appropriate. This is because the bandwidth of Schwinger 
bosons is much smaller than that of the fermion (D <C B) 
and the effects of the quantum fluctuation of bosons can 
be negligible for T>0.5T C where Dt/(2S R T) < 1. The 
memory function is evaluated to leading order in 1/S R , 
and thus, the resistivity is written as 
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where T(Ri) represents the spatial spin correlation de- 
fined as 
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FIG. 2. Calculated resistivity as a function of T/T c for 
several doping concentrations. The results for x = 0.2 and 
x = 0.25 are almost the same and are overlapping. The corre- 
sponding inverse relaxation time (1/t) is shown in the inset. 
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Here, f(ep) is the Fermi distribution function of the spin- 
less fermion, K=pi —j>i is the momentum transfer of the 
fermion due to scattering, and the average (. . .) is eval- 
uated by the mean field Hamiltonian in Eq. (0). The 
same expression given by spin- variables has been derived 
in Ref. [l^. It should be noted that in this model the 
actual value of the resistivity in units of ha/ e 2 does not 
depend on that of t. This is because ep and Kf are scaled 
by t, and r is scaled by 1/t; thus, KfT becomes indepen- 
dent of t. 

Results of the calculated resistivity (p) are shown in 
Fig. |^ as a function of the temperature (T) for several 
doping concentrations. The resistivity monotonically in- 
creases with increasing T in a ferromagnetic state for all 
doping concentrations. In a paramagnetic state, the re- 
sistivity still weakly increases, i.e., metallic for x < 0.2. 
In the case of x = 0.2 to 0.25, the T-dependence almost 
vanishes for T>1.2T C . For a further high doping concen- 
tration, x > 0.2, the resistivity comes to weakly increase 
again. On the other hand, the corresponding inverse re- 
laxation time (1/t) shown in the inset of the figure shows 
a different T-dependence: it decreases in the paramag- 
netic state for all cases of x. This difference clearly in- 
dicates the importance of the T-dependence of Kf since 
p oc l/(KfT). 

These behaviors in the paramagnetic state can he. un- 
derstood in terms of Fisher and Langer's schemejlil al- 
though it was originally proposed to analyze the electron 
transport in the transition metal ferromagnets. The in- 
verse relaxation time in the DE model can be rewritten 
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FIG. 3. An approximate temperature dependence of T(0), 
1/t, and p obtained by Eqs. (10)-(12), respectively, in the 
paramagnetic state for x = 0.3. Each result is normalized by 
the value at T c . 



where f(Ri) is the decaying oscillatory function, and 
D(Ep) is the density of states at the Fermi level. In 
the case of the transition metal ferromagnets, T(0) does 
not depend on T in the paramagnetic state; thus, the 
temperature dependent part of thc_resistivity is wholly 
determined by T{R Z ) for R z ^ OJB In the case of the 
DE model, however, the term with Ri = depends on T 
through C as 



r(0) ~ 8S 2 R {7C 2 
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for S B , > 1 



(10) 



where C is the correlation function between nearest 
neighboring spins. It gives a dominant contribution in 
1/r. Further, D(E F ) (ocl/W/) and Kf also depend on 
T; therefore, the temperature dependence of 1/t and p 
are given by 
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These quantities for x = 0.3 are plotted in Fig. |^ and 
roughly agree with the calculated results shown in Fig. 
||. The rounding of the calculated p at very close to T c , 
however, must come from T(Ri) at Ri^Q. 

It is worth to note that the calculated results in the 
temperature dependence of the resistivity are different 
from those in Ref. [l^, although the memory function for- 
malism is adopted in both cases. The discrepancy is at- 
tributed to the fact that in the present calculation the 
temperature dependence of the electronic structure, that 
is, Wf and Kf (shown in Fig.Jl]), are taken into account, 

as well as that of T,i T {Ri)f{ R i)- In fact, E t r (^)/(-Ri) 
plotted in Fig. |] shows a peak at T < T c which is smeared 
out in the resistivity due to the variation of the electronic 
structure. The peak structure originates from the pro- 
cess in which bosons in the condensate part are scattered 
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FIG. 4. The temperature dependence of ^ T(Ri)f(Ri) for 
several doping concentrations. Results are normalized by the 
value at T c . 

to the non-condensate part, and vice versa.0 It should 
be noted that, since the number of condensate bosons 
is macroscopically large, the lowest order perturbational 
treatment for such scattering processes might overesti- 
mate the scattering amplitude, and the resistivity might 
become slightly smaller than the present results because 
of higher order perturbations. In any case, the resistivity 
is unlikely to show the peak in the ferromagnetic state. 

On the other hand, our results are similar to the pre- 
vious results by Kubo and Ohatafl except for the dis- 
continuity in dp/dT at T c . In Ref. ||, only the shortest 
correlation of scattering (Ri = 0) is included and this is 
a suitable approximation as discussed above. It should 
be noted, however, that the physical mechanism lead- 
ing to the weak temperature dependence of the resis- 
tivity in the paramagnetic state is completely different. 
The mechanism is a result of the cancellation between 
the band- width (W/), the kinetic energy (Kf), and the 
nearest-neighbor spin correlation (C) in a somewhat com- 
plicated manner (Eq. Jl^)) in our calculations. On the 
other hand, these quantities used in Ref. || are temper- 
ature independent in the paramagnetic state. Further, 
the discontinuity in dp/dT at T c comes from the sudden 
change of Wf, which is absent in our results since the 
contribution from T(Ri) with Ri ^ is not neglected 
near T c and Wf is a smooth function of T. Although 
our results are also similar to the results by Furukawa,0 
it is likely to be just a coincidence because the spatial 
correlation of core spins between different sites are ne- 
glected and scattering processes responsible for the re- 
sistivity are different. It should be noted further that 
our results qualitatively agree with the recent results by 
Monte Carlo simulations £3 For a = 4k, at T c , we obtain 
p ~ 2 x 10~ 3 ilem, whose order also agrees with these 
results. 

To conclude, using the Schwinger boson approach, we 
have calculated the resistivity in the double-exchange 



model. In this approach, the fermion bandwidth has been 
determined by the absolute value of the hopping ampli- 
tude giving a physically reasonable temperature depen- 
dence in contrast to the conventional Schwinger boson 
approach. The resistivity monotonically increases with 
increasing temperature in the ferromagnetic state, which 
is different from the previous results in Ref. [l^. In the 
paramagnetic state, the resistivity is dominated by the 
short range correlation of scattering. Although the be- 
havior slightly changes depending on the doping concen- 
tration, the temperature dependence almost vanishes due 
to a cancellation between the nearest-neighbor magnetic 
correlation, the fermion bandwidth, and the fermion ki- 
netic energy. These results imply the importance of the 
temperature dependence of the electronic structure of the 
conduction electron in the both ferromagnetic and para- 
magnetic phases. 

The results agree with the experiments for "higher" 
doping concentrations £~0.3, where experimentally ob- 
served resistivity is relatively low as a whole (p~ 5 x 10 -3 
ficm at T c ) and shows a metallic behavior in the param- 
agnetic state. For lower concentrations x~0.2, however, 
the experimentally observed singular behavior around 
the transition temperature and the insulating behavior 
in the paramagnetic state are not reproduced in our cal- 
culations. Other effects might play a crucial role in co- 
operation with the double exchange mechanism for the 
system with such the lower concentration. 

One of the authors (Ishizaka) would like to thank T. 
Hiroshima for helpful discussions. 
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